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Multigap Superconductivity in Y 2 C 3 : A 13 C-NMR Study 
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We report on the superconducting (SC) properties of Y2C3 with a relatively high transition 
temperature T c = 15.7 K investigated by 13 C nuclear- magnetic-resonance (NMR) measure- 
ments under a magnetic field. The 13 C Knight shift has revealed a significant decrease below 
T c , suggesting a spin-singlet superconductivity. From an analysis of the temperature depen- 
dence of the nuclear spin-lattice relaxation rate 1/Ti in the SC state, Y2C3 is demonstrated to 
be a multigap superconductor that exhibits a large gap 2A//cbT , c = 5 at the main band and 
a small gap 2A/ksT c — 2 at other bands. These results have revealed that Y2C3 is a unique 
multigap s-wave superconductor similar to MgB2. 
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The discovery of superconductivity in MgB2, which 
exhibits a high superconducting (SC) transition temper- 
ature T c ~ 40 K, has attracted much interest. 1 Moti- 
vated by this discovery, intensive effort is devoted to 
the search for a new high-T c material in a similar sys- 
tem that contains light elements B and C. Meanwhile, 
Amano et al. reported that Y2C3 prepared under high 
pressure (~ 5 GPa) is a superconductor with a rela- 
tively high T c ~ 18 K, 2 although the superconductiv- 
ity in this compound was already reported to emerge at 
T c ~ 6—11 K. 3 As for SC characteristics, the specific heat 
measurements on the newly synthesized high-purity sam- 
ples of Y2C3 have revealed that a gap size of the respec- 
tive samples with T c = 11.6, 13.6, and 15.2 K increases 
as 2A//sbX'c=3.6, 4.1, and 4. 4. 4 This result raises a ques- 
tion why T c and 2A//cbX , c vary significantly depending 
on sintering conditions. Further systematic experiments 
are required to gain deep insight into the SC character- 
istics of this compound. 

Y2C3 crystallizes in the cubic Pu2C3-type structure 
(space group /43c?) without an inversion center, consist- 
ing of the dimers of carbon atoms. The SC properties 
of the sample previously reported by Amano et al. did 
not exhibit a single SC transition as seen in the inset 
of Fig. 1(a), pointing to a contamination of extrinsic 
multiple phases. 2 Recently, Akutagawa et al. have suc- 
ceeded in preparing a single phase of Y2C3, which enables 
us to extract intrinsic electronic and SC properties in 
Y2C3. From the specific-heat measurements of this sam- 
ple, the Sommerfeld coefficient 7 ~ 6.3mJ/mol-K 2 and 
Debye temperature 6*d ~ 530 K were estimated for the 
sample with T c = 15.2 K. 4 This result suggests that its 
high Debye temperature makes T c relatively high despite 
its small Sommerfeld coefficient. As in MgB2, the light- 
element constituent like boron and carbon plays a vital 
role in enhancing T c in general. In the SC state, the tem- 
perature (T) dependence of the specific heat exhibits an 
exponential decrease with 2A//cbTc = 4.4 upon cooling 
well below T c , suggesting a strong-coupling isotropic su- 
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perconductivity. From an other context, it is noteworthy 
that a novel SC nature for CePtsSi and Li 2 Pt 3 B without 
inversion symmetry is a recent interesting topic because 
the admixture of spin-singlet and spin-triplet SC state is 
shown to emerge due to the spin-orbit coupling. 5,6 Like- 
wise, determining the order-parameter symmetry and a 
detailed gap structure is an underlying issue in the newly 
synthesized high-quality Y2C3 without inversion symme- 
try. 

In this letter, we report on the SC order-parameter 
symmetry and gap structure of Y2C3 with a relatively 
high T c = 15.7 K (H = 0) via 13 C nuclear-magnetic- 
resonance (NMR) measurements under a magnetic field. 
Y2C3 was synthesized by arc melting and high pressure. 4 
The sample was confirmed to nearly consist of a sin- 
gle phase by X-ray diffraction analyses, with the forma- 
tion of a primitive Pu2C3-type structure. The polycrys- 
talline sample for 13 C-NMR measurement was slightly 
enriched with 13 C ( 12 C : 13 C = 9 : 1) in order to improve 
the NMR signal-to-noise ratio. T c = 15.7 and 12.2 K 
were determined by ac-susceptibility measurements at 
H = and 9.85 T, respectively. The NMR experiment 
was performed by the conventional spin-echo method at 
H = 9.85 T in the T range of 1.8 - 70 K. 

Figures 1(a) and 1(b) show the 13 C- NMR spectra of 
Y2C3 for the previous sample reported in ref. 2 and 
the present sample, respectively, in the normal state at 
T = 15 K and H = 9.85 T. Note that the spectra for the 
previous sample are composed of 13 C- NMR signals aris- 
ing from Y 3 C, YC2, and Y2C3, demonstrating the con- 
tamination of extrinsic multiphases, whereas the spectra 
for the present sample consist of a nearly single peak from 
Y2C3 with a small contamination of YC2. The NMR in- 
tensity for each phase coincides with the x-ray intensity 
as expected. The full width at half maximum (FWHM) 
in the 13 C-NMR spectrum of Y2C3 is as small as 8 kHz, 
ensuring good sample quality. Actually, a single SC tran- 
sition at T c = 15.7 K was corroborated by the suscepti- 
bility measurement at H = 10 Oe as seen in the inset of 
Fig. 1(b). 
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Fig. 1. (Color online) 13 C-NMR spectra of Y2C3 for (a) the pre- 
vious sample reported in ref. 2 and (b) the present sample in the 
normal state at T = 15 K and H = 9.85 T. Note that the spectra 
for the previous sample are composed of NMR signals arising 
from Y3C, YC2, and Y2C3, demonstrating the contamination of 
extrinsic multiphases, whereas the spectra for the present sample 
nearly consist of a single peak from Y2C3 with a small contami- 
nation of YC2 • The insets of both figures show the T dependence 
of SC diamagnetic susceptibility down to 4.2 K. 



Figure 2(a) shows the T dependence of 13 C Knight 
shift (KS) in Y2C3, which is determined relative to the 
resonance frequency of tetramethylsilane (TMS) as a ref- 
erence substance (iC[TMS] ~ Oppm). A clear decrease in 
KS and an increase in FWHM below T c are indicated in 
Figs. 2(a) and 2(b), which are derived from the T depen- 
dence of NMR spectra as shown in the inset of Fig. 2. 
In intermediate fields (H c i <C H <C H C 2) where the vor- 
tices form a dense lattice, we estimate coherence length 
and the distance between the vortices to be d ~ 160 A, 
and £ ~ 34 A, respectively. 8 As a result, a diamagnetic 
field is led to be Tfdia ~ — O.3 0e using the relation 
H dm = -i/ c iln(/3e- 1 / 2 d/0/lnK, 7 ' 8 (k = A/0 and hence 
a diamagnetic shift is obtained as ~3x 10 -4 % at 
ff=9.85T. Here, we used H cl = 3.3 mT, 4 the London 
penetration depth A = 4470 A, 4 /3 — 0.381 for the case of 
triangular lattice, 7 d ~ 160 A, and £ ~ 34 A. Thus, the 
estimated value of K^ a is one order of magnitude smaller 
than the decrease in KS observed below T c , demonstrat- 
ing that the decrease in KS is due to the reduction of 
spin susceptibility associated with the onset of the spin- 
singlet SC state. If the spin susceptibility was assumed 
to vanish at low T due to the formation of spin-singlet 
Cooper pairing, the orbital and spin part of KS are tenta- 
tively estimated to be K oth ~ 0.028 % and K s ~ 0.005 %, 
respectively. A possible cause for the increase in FWHM 
may be due to an inhomogeneous distribution of vortex 
lattices, which eventually makes either d or A distribute. 
Further systematic NMR measurements at low H are re- 
quired for inspecting a structure of vortex lattices. 

Next, we deal with the T dependence of nuclear spin- 
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Fig. 2. (Color online) T dependences of (a) 13 C Knight shift and 
(b) full width at half maximum (FWHM) of 13 C-NMR spectrum 
in Y 2 C 3 at H = 9.85 T. T c = 12.2 K is determined by the ac- 
susceptibility measurement at H = 9.85 T as shown by an arrow 
pointing downward. The inset shows the NMR spectra at T = 
2.9, 5.4, 7.8, 9.5, 12, and 18 K, respectively. 



lattice relaxation rate (1/Ti) in order to clarify the gap 
structure. The 1/Ti for 13 C with nuclear spin 7 = 1/2 
is uniquely determined from a simple exponential recov- 
ery curve of nuclear magnetization given by the relation 
[M(oo) - M(t)]/M(ao) = exp(-t/Ti). Here, M(t) and 
M(oo) are the nuclear magnetizations at a time t after 
the saturation pulse and at the thermal equilibrium con- 
dition, respectively. Figure 3 presents the T dependence 
of 1/Ti at H = 9.85 T. In the normal state, the law 
TiT = const, is valid down to T c . In the SC state, X/T\ is 
also precisely measured from a simple exponential recov- 
ery curve such as [M (oo) - M[t)]/M{oo) = exp(-t/T 1 ), 
as seen in the inset of Fig. 4(a). This is because the possi- 
ble contribution to 1 /T\ arising from normal vortex cores 
is very small, if any, when £ - 34 A<< d ~ 160 A. 

The inset in Fig. 3 shows (TiT) const ./(TiT) vs T/T c 
for Y 2 C 3 at H = 9.85 T and for MgB 2 with T c = 29 K 
at H = 4.4 T. 10 Here, (TiT) const . denotes constant val- 
ues in normal state. In the SC state, we note that a 
tiny coherence peak is observed in 1/Ti just below T c for 
Y2C3 as in MgB2. This is indicative of a full gap open- 
ing in Y2C3 as in MgB2. A reason why the coherence 
peak is depressed in these compounds may be due to 
a strong electron-phonon coupling that causes the large 
life time broadening of quasiparticles induced by ther- 
mally excited phonons as reported in ref. 10. In fact, the 
strong-coupling BCS superconductor such as TlMogSey.s 
does not show a clear coherence peak. 11 Note that the 
T dependence of 1/Ti well below T c does not exhibit a 
simple exponential decrease, but seems to have a kink at 
around T = 5 K. In order to gain further insight into this 
unique and relevant relaxation behavior with a possible 
gap structure in the SC state for Y2C3, we present in 
Fig. 4 the Arrhenius plot of (TiT)/(TiT) const . vs T c /T 
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Fig. 3. (Color online) T dependence of 1/Tj for Y2C3 with 
T c = 12.2 K at H ~ 9.85 T. A tiny coherence peak just be- 
low T c is observed for Y2C3 as in MgB2 as shown in the inset. 
The inset shows the (TiT) const ./(TiT) vs T/T c curve for Y2C3 
(solid circles) at H = 9.85 T and for MgB2 (solid squares) with 
T c = 29K at H = 4.4T. 10 Here, (TiT) const . denotes constant 
values in normal state. 

with T c = 12.2 K at i? = 9.85 T. It is evident that a 
power-law behavior in 1/Ti such as \fT\T oc T 2 (see 
the dashed line in the figure) is not valid at all. Instead, 
when noting that a line in this plot corresponds to an 
exponential T dependence in 1/T\T, it is supposed that 
a large full gap seemingly opens in a high-temperature 
regime in the SC state, but low-lying quasiparticle ex- 
citations in a low-temperature regime are dominated by 
the presence of a small full gap. In fact, from respective 
slopes in this plot in the T range of 5K-T C = 12.2 K and 
at temperatures lower than 5K, gap sizes are estimated 
to be 2A/ksT c ~ 5 and 2, which suggests that two kinds 
of SC energy gaps exist, namely, multigap superconduc- 
tivity takes place in Y2C3. We stress that this novel re- 
laxation behavior in the SC state for Y2C3 is not due 
to some inhomogeneous effect originating from the pres- 
ence of vortex cores and/or a distribution of T c because 
1 /T\ is uniquely determined from the simple exponential 
recovery curve of nuclear magnetization as shown in the 
inset of Fig. 4. 

Here, we apply a phenomcnological multigap model 
for nodeless superconductivity to understand the novel 
relaxation behavior in the SC state. Figure 5 shows the 
T dependence of \jT{F. In such a model, Ti(T c )/Ti(T) 
is expressed as 

Ti(T c ) _ a 2 1_ f? J_ 

Ti(T) a 2 + /3 2 T 1 y a) a 2 +/3 2 T X K 0h 

where a and p are defined as the respective fractions of 
N(E F ) x A M with SC gap A a and and a + (3 = 1. 
7V(Ep) and Am are the density of states (DOS) at the 



Fermi level and the hyperfine coupling constant, respec- 
tively. Here, 
1 2 f°° 

_(A) = — - / dE[N 2 (E) + M 2 (E)]f(E)[l-f(E)}, 
Ti k B T c J a 

where N S (E) is the DOS, M S (E) is the anomalous DOS 
originating from the coherence effect inherent to a spin- 
singlet SC state and f{E) is the Fermi distribution func- 
tion. 12 Note that N S (E) and M S (E) are averaged over an 
energy broadening function assuming a rectangle shape 
with a width 26 and a height 1/2S. 13 We use 5/A(0) = 
0.3 in the calculation. A theoretical curve based on the 
multigap model is actually in good agreement with the 
experiment using j3 = 0.75 and 2Aj3/k B T c — 5 for 
the main band, and a = 0.25 and 2A a /k-QT c = 2 for 
other bands as shown by the solid curve in Fig. 5. It 
is notable that the large gap at the dominant Fermi 
surface is larger than the weak-coupling BCS value of 
2A/k B T c = 3.5, indicating a strong electron-phonon cou- 
pling and being consistent with the specific- heat result. 
The present 13 C-NMR has revealed that the supercon- 
ductivity in Y2C3 is characterized by a large gap at the 
main Fermi surface and a small gap at others. This may 
be consistent with the band calculation which shows the 
presence of Fermi surfaces consisting of three dimensional 
multishects due to the hybridization between Y-d de- 
rived states and antibonding C-dimers derived p-states. 9 
We should pay attention to the relationship between the 
multigap and T c . Although Y2C3 is a superconductor 
with no inversion symmetry, the present experiments 
have revealed that this noncentrosymmetric compound 
is a spin-singlet superconductor with full gaps at all the 
Fermi surfaces, and hence rules out the possibility of the 
admixture of spin-triplet order parameter which is the 
recent underlying topic for the superconductors with no 
inversion symmetry. 




T /T 

Fig. 4. (Color online) Arrhenius plot of (TiT)/(TiT) const, vs 
T c /T with T c = 12.2 K at H = 9.85 T. From the respective slopes 
in the T range 5 K — T c = 12.2 K and at temperatures lower than 
5K, gap sizes are estimated as 2A/k^T c ~ 5 and 2. The inset 
shows a simple exponential recovery curve of nuclear magnetiza- 
tion given by the relation [M(oo) - M(t)]/M(oo) = exp(-t/Ti). 
Here, M(t) and M (oc) are the nuclear magnetizations at a time 
t after the saturation pulse and at the thermal equilibrium con- 
dition, respectively. Even in the SC state at T = 7K, note that 
1/Ti for 13 C with nuclear spin / = 1/2 is uniquely determined 
(see the text). 
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Fig. 5. (Color online) T dependence of 1/TjX for Y 2 C 3 with 
T c = 12.2 K at H = 9.85 T. A phenomenological multigap 
model for nodeless superconductivity is applied to understand 
the novel relaxation behavior in the SC state. The solid curve is 
a theoretical curve based on the multigap model with (3 = 0.75 
and 2A / 3/A;bT c = 5 for the main band, and a = 0.25 and 
2A a /kftTc = 2 for other band (see the text). The inset shows 
the T dependences of the order parameters A a and Ag. 



In conclusion, the superconducting properties of Y2C3 
with a relatively high transition temperature T c = 15.7 K 
have been investigated using the 13 C nuclear- magnetic- 
resonance (NMR) method under a magnetic field. The 
Knight shift has revealed a significant decrease below T c , 
suggesting the spin-singlet superconductivity. The nu- 
clear spin-lattice relaxation study in the SC state has 
revealed that Y2C3 is a multigap superconductor that 



exhibits a large gap 2A/fceT c = 5 at the main band and 
a small gap 2A/k^T c = 2 at others. These results have 
revealed that Y2C3 is a unique multigap s-wave super- 
conductor similar to MgB2. 
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